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LA W OF ERROR IN THE POSITION OF A POINT IN SPACE. 



BY E. L. DE FOBEST, WATEBTOWN, CONN. 

The law of facility of error in the observed position of a point varies 
according as the point is supposed to lie on a given straight line, or in a 
given plane, or in unrestricted space of three dimensions. Most quantities 
which are the subjects of observation are capable of having their magnitude 
represented by the length of a straight line, and their errors are represented 
by the errors in the position of the point which forms one extremity of the 
line, the other extremity being fixed. Taking this line as an axis of X, the 
probability y that the point of error will fall within any arbitrary but very 
small interval dx on this line, is a function of the form 

y = *%-«-»", (1) 

J/3T 

where x is the distance from the most probable position of the point to the 
middle of the interval dx, and A is a constant. The most probable position 
of the point, and the value of h, are to be determined by known methods. 
If an observed point may be anywhere on a given plane, as in the case of 
shot-marks on a target, the probability that the point will fall within any 
very small area dxdy whose middle point has the coordinates x and y, will be 

z ^1^2 ( ^ < fy e -Wfr'-HW') j (2) 

The positions of the axes of X and Y, and the values of h x and h 2 , are 
to be found as in my recent article, Analyst, Jan. to May '81. In this 
and other articles (Sept. '81, p. 145), it was shown that the functions (1) 
and (2) are the limiting forms of the series of coefficients in the expansion, 
to an infinitely high power, of a polynomial of one and two variables 
respectively. On similar principles, the law of errors in space of three 
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dimensions may be obtained by the expansion of a polynomial of three va- 
riables, whose coefficients are all positive and their sum equal to unity. It 
is the object of the present paper to show this, and to demonstrate the lim- 
iting form of the expansion. 

Taking at first the true position of the observed point as an origin, and 
drawing through it any three rectangular axes, we will refer to these axes 
any other positions which the point may take through errors of observation. 
The coordinates x, y, z of any such point of error will be regarded as mul- 
tiples of the arbitrary intervals or units of measure Ax, Ay and Az respect- 
ively, which may be taken as small as need be. If for example a single 
observation places the point at x = 4Ax, y = 2Ay, z = SAz, we will write 
L 4 2 , 3 for the probability that this particular error will occur, and in gen- 
eral, we write L a-b e for the probability of occurrence of the error x — aAx, 
y—bAy, z = cAz, where a, b, c are any whole numbers, positive or negative. 
Now take the probabilities L as coefficients, and their sub-indices as expo- 
nents off, 7], £, in the polynomial 

Under this notation, the terms of the polynomial are supposed to be ar- 
ranged in the form of a rectangular block, so that exponents and sub-indices 
vary from term to term by differences equal to unity. Each term is of the 
form included within the brackets, and the signs of summation denote that 
the whole polynomial is the sum of all the terms, and is formed by assigning 
to one of the exponents and sub-indices, a for instance, all the integer val- 
ues from — m to m, while b and c remairi constant, then in this sum assign- 
ing integer values to another exponent, b for instance, from — m to m, while 
c remains constant, and finally in this sum assigning integer values to c 
from — m to m. In other words, the block-formed polynomial is the sum 
of all the different terms which can be formed by giving to a, b and c, in- 
dependently of each other, all the integer values from — m to m. We might 
have adopted a notation analogous to that which I employed for polynom- 
ials of two variables, in Analyst, Jan. '81, p. 4, but it would be quite 
cumbrous, if the eight corner terms of the block were written in full. 

The number m is supposed to be taken so large that no point of error 
will occur whose coordinates exceed the limits x = ± mAx, y = ± mAy, 
z = ± mAz. Any entire polynomial of three variables can be put in this 
block form, by adding or intercalating terms with zero coefficients if req'd. 
The whole number of terms in the block is (2m -f l) 3 . From the manner 
in which the coefficients in the successive powers of a polynomial are formed, 
it is evident that when n observations are made, the probability that the 
sum of the errors in the x direction will be sAx, and that at the same time 
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the sums of those in the y and z directions will be tAy and v Az respectively, 
is the coefficient of £ 'rf£° in the expansion of the polynomial (3) to the nth 
power. 

If L is the probability that a particular error or event will occur, in one 
observation, then in n observations the most probable number of times for it 
to occur is nL. (Analyst, May '79, p. 66.) This is strictly true when 
nL is a whole number, and approximately true in other cases. Hence, the 
most probable total effect, in the x, y and z directions, of the possible error 
4Jx, 2Jy, Zdz for example, will be 

nL i, 2, z( iJx )> wi 4, 2, s( 2 4/)> nL i,2,s( 8Jz )> 

and in like manner for the other errors. The most probable algebraic sum 
of the n errors which occur will be, in the x direction, 

nA^ t =^2i=%,2%=^JflL a . lit )\ (4) 

where the coefficient a takes the same values as the sub-index a does. This 
is the approximate value of the exponent of c in that term of the expansion 
of (3) whose coefficient is a maximum. The most probable value of the 
arithmetical mean of the n errors which occur, in the x direction, is found by 
dividing (4) by n. The quotient thus obtained is evidently the statical 
moment, with respect to the YZ plane, of the system ot coefficients L re- 
garded as the masses of material points. It also represents the lever arm 
of the system with respect to that plane, since 1L = 1. In like manner 
it will be found that the most probable arith. means of the errors in the Y 
and Z directions, respectively, are the lever arms of the system with respect 
to the XZ and XY planes. Thus the most probable mean position of the 
n points of error is the centre of gravity of the system of masses L. 
From what was shown in my article of May '80, p. 81, respecting the po- 
sition of the centre of gravity of the coefficients in an expanded polynomial, 
joined with the fact that the coefficients in the expansion are not altered by 
subtracting a constant quantity from each exponent in the first power, so 
as to bring the place of f ^ (or x^z" in the notation there used) into the 
middle of the block instead of at one corner, it follows that the lever arms 
of the system of coefficients in the expansion of (3) to the nth power are n 
times what they are in the first power, the coordinate planes remaining the 
same. Hence, according to (4) <*fcc , the maximum coefficient in the expan- 
sion coincides approximately with the centre of gravity of the whole sys- 
tem of coefficients in the expansion. It will otherwise appear hereafter, 
that this is true at the limit, when n becomes infinite. 

Let us write the expansion of (3) to the nth power as follows, retain- 
ing the same block form. 

77 Vc = mn yb~mn Va = jnn // ta-JtycX /n 

u •'l=-« i l=-i»^«=-ni^,i 1 |t fj C )> \°) 
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where a, b, c take only integer values as before. The whole number of 
terms in the expansion is (2mn-fl) 3 . We shall now employ the method of 
indeterminate coefficients, to find a relation between the (2m + l) 3 coeffi- 
cients L in the first power, and any similar block of an equal number of 
coefficients I in the expansion. Since the L and I are independent of f, rj 
and f , we may take the relation 

u" = U, (6) 

and differentiate it with respect to the variables, obtaining 

which being divided by (6), gives 

•"©-(if)- -OHf)' » p (|)=«(f)- P) 

In the first of these three equations (7) substitute the expressions for u 
and Z7as in (3) and (5), and also expressions for -^ and — ^- as obtained 

from the same by differentiation. This gives an equation, each member of 
which contains the product of two block-formed polynomials. Any such 
product will itself be a block-formed polynomial, and any term in it will 
have for its coefficient the sum of products of L and I, such that this sum 
can be expressed in the block notation. Forming in this way the coefficient 
of £ i ~ 1 r/£'° in the product in each member of the equ'n referred to, and plac- 
ing them equal to each other by the principle of indeterminate coefficients, 
we get 

l[^t = ^> = -«^« = -'»( a J-'-a,-b,~c '(i+«),(j+*),(H«)lJ 

=^5=-)»^i = -m^o=-«lL(* - r" a )J-'-a,-b -c h'+a)-(j+ b),(k+c) I - (°) 

In the second member of this equation, let the part which does not have 
the coefficient i be transferred to the first member, then change the signs of 
both members and divide by 71 + 1. Take Fas an auxiliary letter and write 

' — ^l = -vi^-b = -m^'a = ~m(J-'-a,~b,-c'(i+a),U+l'), (*+<=))• (.") 

Then (8) becomes 

c = m bznm av=.m i , 1 

I 2' 2 (ai_ a ,_ 6 ,~^(i+a), (j+S),(ft+c)) = — j-- T V. (10) 

In like manner, starting from the two last equations in (7), and equating 
to each other the coefficients of f^-if* in the one case and those of £*r/£ k ~ l 
in the other, we get after reduction 

czzin, h~m a~m / \ —7* 

c z= — m 6 — — «i a = — *» \ W "T* -*■ 

2 2" 2" (ci_ ,_5,_c W).a+»), (*+«)) == rxT 17 * ( 12 ^ 

•= _mS = — mo = — m\ / 7* -J- i. 
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It will be seen that the terms in (10), (11) and (12) all contain products 
of L and I, so taken that any L occupying a given position in its block, is 
the multiplier of that I which occupies just the opposite position in the equal 
block of terms in the expansion. These three equations express the rela- 
tion which exists between the coefficients L in the given polynomial, and 
any similar block of coefficients I in the expansion to the nth power. The 
middle coefficient in this block is lij, k , and its coordinates are iAx,jAy and 
kAz, the expansion being referred to the same axes as the first power. The 
coefficients I are also regarded as the masses of material points. When the 
exponent n is made very large, or infinite, the masses I attain their limiting 
values, and we denote by w the limiting value oil iik . Supposing the p'ts 
to be set so close together as to become consecutive, Ax, Ay, Az are repre- 
sented by dx, dy and dz, the coordinates of the point whose mass is w are 

x = idx, y — jdy, z = kdz, (13) 

and we wish to express w as a function of x, y and s. Let n be made an 
infinity of the second order. The system of material p'ts will then extend 
throughout infinite space, their number (2mn. + l) 3 being the whole number 
of consecutive points in the infinite block, while (2m -f- l) 3 is the number 
included within the small block under consideration which has w at its cen- 
tre, so that the portion of space occupied by this last block is infinitesimal 
in its length, breadth and depth. The masses of the consecutive points 
within its limits may therefore be regarded as varying from each other by 
constant differences, under the assumption, of course, that w is a continuous 
function of x, y, z. 

Let djc, d v w and d z w denote the differentials of w taken in the three co- 
ordinate directions. The values of all the (2m + l) 3 masses I in the block 
will be formed from w by successive additions and subtractions of its differ- 
entials, and will be the terms in this block-formed total, 

2tZ™ m 2l=™ m Z:='l m (w+ad x w+bd y w+cd z w). (14) 

Substitute them for the correspond'g terms I in (10), (11) and (12). Collect 
separately the coefficients of w, d x w, d y w and d x w in the result, remembering 
that IL = 1. Employ auxiliary letters, such that a v u 2 , a 3 denote the 
three sums of all the products formed by multiplying each L by its first, 
second, and third sub-index respectively; also let fl lt /9 2 , /3 S denote the 
three sums of the products of each L into the square of its first, second, and 
third sub- index respectively; and let y lt f 2 , y z denote the three sums of 
the products formed by multiplying each L by the product of its first and 
second, first and third, and second and third sub-indices respectively. It 
will be found that (10), (11) and (12) are then reducible to 
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— a x w-\- ^ i d x w-\-y 1 d t/ w-\-y 2 c] z w = Aw— a x d x w — a 2 d y w — a 3 d z iv), J 

n-f-1 I 

—azw+r^djvl P 2 d v w + r 3 djD — —^(w-a 1 d x w—a 2 d y w-a 3 d z to) ) )■ (15) 

— k I 

— a z w+r 2 d x w+r z d y w+p z d,w = -^-r{v)—a 1 d x w—a 2 d y w — a z d.w). 

These are differential equations in which a, , /9j &c. are constants. From 
the manner in which the latter are formed, it is evident that if the coeffi- 
cients L of the given polynomial are regarded as the masses of material 
points, the intervals between them being Ax, Ay, Az, represented by dx, dy, 
dz at the limit, then a 1 dx, a 2 dy, a z dz are the statical moments of the sys- 
tem with respect to the coordinate planes YZ, X£and XY. Also /? 1 (c/a;) 2 , 
fi 2 [dyf and ft z {dzf are the three sums of the products formed by multiply- 
ing each L into the square of its distance from the three planes respectively. 

Likewise y x dxdy, y 2 dxdz and y z dydz are the three sums of the products 
formed by multiplying each L into the product of its distances from the 
planes FZand XZ, YZ a,nd XY, and XZ audi XY respectively. It is evi- 
dent that if the coefficients L in the given polynomial (3) are such that their 
centre of gravity falls at the origin, that is, at the place of i , , , the stati- 
cal moment of the system with respect to any plane through this centre will 
be null, and we shall have a t = 0, a 2 = 0, a 3 =0. Furthermore, if the 
coefficients and assumed coordinate axes are such that these axes are the free 
axes of the system of masses L, that is, the principal axes through the cen- 
tre of gravity, we shall have also j' 1 =0, y 2 — 0, 7" 8 =0. This follows from 
the known properties of the principal axes in bodies of three dimensions. 
(See for instance Chapter II. of Vol. II. of Poisson's Traite de Me'oanique.) 

Of course, if the possible points of error and their probabilities are arbi- 
trarily distributed, in other words, if any arbitrary set of values is assigned 
to the coefficients L in (3), their centre of gravity will not in general coin- 
cide with the origin or place of L w , , and the axes taken will not in gen- 
eral be free axes. Bug the position of the free axes can always be found by 
means of the formulas demonstrated in Poisson's treatise, and being at right 
angles to each other, they can be taken as coordinate axes. By reasoning 
precisely analogous to that employed in my article on errors in two dimen- 
sions (Analyst, March '81, pp. 44 to 47), it will appear that when the 
coefficients L are referred to any new rectangular axes, and the coefficients 
I in the expansion are also referred to the same axes, the coefficients in the 
expansion are unchanged and retain the same positions relatively to each 
other which they had under the old system, differing only in position rela- 
atively to the axes. 
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For the sake of simplicity we will assume the same arbitrary interval or 
unit of measure in all directions, so that Ax = Ay = Az, A given system 
of coefficients located in space will be here said to be referred to any given 
rectangular axes, when the variables £ , -q, £ attached to each coefficient re- 
ceive exponents equal to the three coordinates of the coefficient, expressed 
in units of measure. And coefficients whose attached variables have any 
given exponents are said to be referred to given axes, when they are so 
placed in space that their coordinates, expressed in the units of measure, are 
equal to the corresponding exponents. 

Suppose that the directions of the axes in (3) are changed while the ori- 
gin remains the same. Let x lt y if z x and x 2 , y 2 , z 2 be the coordinates, 
in units of measure, of two given points of error referred to the old axes, 
and let x' lt y\, z\ and x' 2 , y' 2 , z' 2 be those of the same points referred to 
the new axes. Multiplying together the two corresponding terms of the 
polynomial, omitting their coefficients L since the products of these have 
the same value in both cases, we find in the two systems 

p'lrrtMXfy'irt _ £*'i+*'jjjff'i+tf's^»'i-H' s> j ^ 

The exponents of the second members are the coordinates of the place of 
the coefficient of the product of the two terms, in the two systems. The 
known relations between the two systems of coordinates give 

aj x = a x\+b y\-\ro z\, x 2 = a x' 2 +b y\+c z' 2) ~| 
y x = a'x' 1 +b'y' 1 +e'z' 1 , y 2 = a'x' 2 +b' y' 2 +c' z' 2 , V (17) 
z x = a"x' 1 +b"y' 1 +c"z' lt z 2 = a"x' 2 +b"y' 2 +c"z' 2> ) 

where a, b, c &c. are the direction-cosines of the new axes, that is to say, a, 
b, e are the cosines of the angles which the new X, Y and Z axes make with 
the old X axis; a', b', c' are those of the angles which they make with the 
Faxis; and so on. By addition we get 

x x +x 2 = a (x\+x' 2 ) + b (y'i+y'J+o (s^ +«',), ) 

2/1 + 3/2 = o'(«'i+« , s) + 6'{y'i+y'a) + o' (a', +*',), } (18) 
z x + z 2 =a"(x' 1 +x' 2 ) + b"(y' 1 +y' 2 )+c'(z' 1 +z' 2 ). ) 
These relations, being of the same kind as those in (17), show that just 
as x lt y x , Zj and x' v y\, z\ are coordinates of the same point in the two 
systems, so the point whose coordinates in the old system are x x -\-x 2 , y, + 
y 2 , z x -f- z 2 , is the same point whose coordinates in the new system are x\ 
-\-x' 2) y\ -\-y' 2 , z\ +z' 2 . In other words, the product keeps the same ab- 
solute position under either system. And as this is true for all the partial 
products which make up the total, it follows that whatever new directions 
we give to the rectangular axes, the origin being unchanged, the absolute 
places of the coefficients in the squared polynomial will be the same. 
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In like manner it may be shown that they will continue the same for 
higher powers, and up to the limit. If we also transfer the origin to some 
new point, and draw axes through it parallel to the new directions as above, 
and denote by a lt b lt o 1 the coordinates of the old origin referred to the 
new, then in the first power of the polynomial all the exponents of f, jy, £ 
will be increased by the constant quantities a lt 6 1 , c x respectively, in the 
second power they will all be increased by 2a lt 2b 1} 2c x respectively, and 
in the nth power, by na 1 , nb^, nc 1 respectively. Thus the positions of the 
coefficients in the expansion, relatively to each other, will be unaffected by 
a change in the origin, or in the directions of the axes, or in both together. 

Supposing now that all the exponents and corresponding sub-indices in 
the polynomial (3) have been altered so as to refer to the new axes, it is ev- 
ident that they will in general be fractional numbers, and »', /, and h likewise, 
whereas we regarded them as integers, in our demonstration of the results (10), 
(11), (12). But we may imagine that the given block of coefficients L in 
(3) is enclosed within a cubical block whose centre is at the desired new or- 
igin and whose sides are parallel to the new axes, and that this new block 
is subdivided into very small cubes by equidistant planes parallel to its 
sides and taken very near together. 

The points at which these planes intersect are supposed to be each occu- 
pied by a coefficient L, only all these coefficients are zero except at the 
points which come nearest to the coefficients iD the original block; and to 
these points the original coefficients are supposed to be transferred. It is 
evident that by diminishing the intervals between the planes, these positions 
of the coefficients L in the new block can be made to approach to their true 
positions in the old block, within any assignable limits. Hence the rela- 
tions (10), (11), (12) will hold good, it being understood that m is now the 
number of new intervals from the centre to the surface of the new block, in 
any one of the three coordinate directions, and that i,j, h denote the num- 
ber of such intervals from the centre of any equal block of coefficients in 
the expansion, to the respective coordinate planes. But we may restore, if 
we please, the larger unit of measure Ax = Ay = Az, by dividing (10), (11) 
and (12) throughout by the ratio which the small interval or unit bears to 
the larger one. The effect of this will be, that the coefficients from m to 
— m, as well as i, j, k, will now, in general, be fractional numbers. The 
same is true of the coefficients from m to — m in (14). Here d x w, d y w and 
d,w are the increments of to corresponding to those increments of x, y and 
2 which are equal to Ax = Ay = Az, and which are represented by dx 
= dy = dz at the limit. 

[To be continued.] 



